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Abstract- In this paper an optimal integral controller is proposed 
for load frequency control problem of multi-area electrical power 
system. The design of optimal controller uses modern optimal 
control theory. A state model of a four -area interconnected 
power system is developed. Dynamic analyses have been done 
without controller, with integral controller and with optimal 
integral controllers using Matlab Simulink- Workspace. A 
comparison between some recent trend controllers and proposed 
optimal integral controller is presented further and it has been 
shown that the performance of the proposed controller is 
superior in terms of overshoot and settling time. 
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I. INTRODUCTION 

Large-scale power systems are normally composed of 
control areas or regions representing coherent groups of 
generators. Various areas are interconnected through tie-lines. 
Tie-lines are utilized for contractual energy exchange and 
provide inter-area support in case of abnormal conditions. 
Real and reactive power generations must change accordingly 
to match the load perturbations. Load frequency control is 
essential for successful operation of power systems to 
maintain frequency within the limited band [1]. To accomplish 
this, it becomes necessary to automatically regulate the 
operations of main steam valves or hydro gates in accordance 
with a suitable control strategy, which in turn controls the real 
power output of electric generators. The controlling of the 
output of electric generators in this way is termed as automatic 
generation control (AGC) [2]. Automatic generation control is 
the regulation of power output of controllable generators 
within a prescribed area in response to change in system 
frequency, tie-line loading, or a relation of these to each other, 
so as to maintain the schedule system frequency and establish 
interchange with other areas within predetermined limits. 
AGC can be sub-divided into primary and secondary control 
modes. The overall performance of AGC in any power system 
depends on the proper design of both primary and secondary 
control loops. Secondary controllers are designed to regulate 
the area control errors to zero effectively [2]. So the overall 
performance of AGC in any power system depends on the 
proper design of both primary control loop (selection of R) 
and secondary control loops (selection of gain for 
supplementary controller). A control strategy is needed that 
not only maintains constancy of frequency and desired tie- 
power flow, but also achieves zero steady state error and 
inadvertent interchange. Among the various types of load 
frequency controllers, the most widely employed are integral 
(I), proportional plus integral (PI) and proportional plus 
integral plus derivative (PID) controllers. The use of modern 
optimal control theory can improve the dynamic performance 
of these controllers [3]. Over the past decade many control 



strategies for load frequency control of power systems such as 
linear feedback [4], variable structure control [5] etc. have 
been proposed in order to improve the transient response. An 
adaptive controller with self-adjusting gain setting was 
proposed for LFC in [6]. Artificial Neural Networks have been 
successfully applied to the LFC problem [7], [8]. Moreover 
fuzzy logic control techniques for LFC problem are mostly 
based on fuzzy gain scheduling of integral [9] and 
proportional integral (PI) [10], [11], [12] controller parameters. 
Recently hybrid neuro fuzzy control strategy is proposed for 
LFC to improve performance of the fuzzy controller [13]. 
Recently many researchers have applied Genetic algorithm [14] 
and neural networks [15] to controllers design to improve the 
dynamic performance of the power system. In a control 
schemes LFC by automatic design of the fuzzy rules in the 
fuzzy gain scheduling control approach by genetic algorithm is 
proposed in [16]. 

This paper uses the optimal control strategy to optimize 
integral controller for four-area interconnected power system. 
Simulation results shows that the proposed controller in terms 
of dynamic performance such as peak overshoot, settling time 
and integral absolute error is better than integral controller and 
many intelligent controllers used recently. 

II. FOUR- AREA POWER SYSTEM 

An interconnected power system is divided into control 
areas connected by tie lines. A four-area interconnected power 
system (Fig.l) of a non-reheat thermal plant is used to explain 
the motivation of the proposed method. The frequency 
deviation in all areas severely affects the production and the 
quality of frequency sensitive industries such as the spinning 
and weaving industry, petrochemical industry, pulp and paper 
industry, etc. Furthermore, the lifetime of machine apparatuses 
on the load side will be reduced. The tie-line power flow and 
frequency of the area are affected by the load changes. 
Therefore, it can be considered that each area needs its system 
frequency and tie-line power flow to be controlled. 
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Fig. 1 Four interconnected control areas (six tie lines) 



III. UNCONTROLLED FOUR-AREA POWER SYSTEM 

A four-area power system connected by six tie lines 
illustrated in Fig. 1 is taken as a test system in the study. Each 
area decides steam turbine contains a governor and a generator. 
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x 18 = -2?r7" 42 x 4 + 2^-7 42 x 1H (18) 
These eighteen equations can be organized in the following 
vector matrix form: 

x = Ax + Fw 

Where, A is system matrix and F is disturbance distribution 
matrix, 

x = [x; . . . x 18 ] T = state vector and 

AP ( = w = [w, w 2 w 3 w 4 ] r =disturbance vector. 

Using equations x l ,x 2 ,...,x lg , matrix A is formulated, as 
given. 
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Disturbance matrix F is given as: 
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Considering Af i, Af 2 , Af 3 and Af 4 (i.e. change in frequency) as 
system's output variables, the output matrix can be given as: 
y = Cx + Dw 
Where, D=0 and C & x are given as: 

1 00000000000000000 
000100000000000000 
000000100000000000 
000000000100000000 

x r = [x 1 □ □ □ x lg ] 



c= 



The Matlab commands [y,x]=step(A,Fw,C,D,l,t) and plot(t,y) 
are used to see the dynamics of Af i, Af 2 , Af 3 and Af 4 of four- 
area LFC controller for change in demand in w, (AP d i). The 
results are verified through Simulink. 

IV. INTEGRAL CONTROLLER FOR FOUR-AREA SYSTEM 

The block diagram of four-area system with integral 
controller is given in Fig. 2. 

Mathematical Model Development 

The following modifications and additions will take place 
in the model developed without controller: 

1 1 1 (19) 
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(23) 
(24) 
(25) 
(26) 

These twenty two equations can also be organized in the 
matrix form as follows: 



"^"20 ^2^4 ^12~^13 ^32~^17 ^42"^"18 

x 21 = & 3 ^ 7 + x 14 + x 16 + x 17 



b A x w + x l5 a 34 x 16 + x ls 



D=0 and C & x are given as: 

"1 00000000000000000000 
0001000000000000000000 
0000001000000000000000 
0000000001000000000000 
0000000000001000000000 
00000000000001000000 
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x = Ax + Bu + Fw 
Where, 

B is the input matrix and u is control vector. 
Further considering the following, 
m, = -k n x l9 = -k n J" A CE t clt = AP C1 (s) 
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The state space model can be represented by the following 
equation: 

x = Ax + Fw (28) 
Where, the system matrix A can be formulated, as given on 
next page. 

Considering Afi, Af 2 , Af 3 and Af 4 (i.e. change in frequency) 
and AP^ 12 , AP tle> 41 (change in tie line power) as system 
output variables, the output matrix can be given as: 
y = Cx + Dw (29) 



The dynamic responses of Afi, Af 2 , Af 3 and Af 4 of four-area 
LFC controller for change in demand in w, (AP d .) are obtained 
through Matlab program and the results are verified through 
Simulink. 

V. OIPTIMAL INTEGRAL CONTROLLER FOR FOUR-AREA 
SYSTEM 

Optimal control is a branch of modern control theory that 
deals with designing controls for dynamic systems by 
minimizing a performance index that depends on all system 
states. The object of optimal control regulator design is to 
determine the optimal control law w, which can transfer the 
system from its initial state to the final state so that a given 
performance index is minimized. 

w = -Kx (30) 
Here K is feedback gain matrix. 

A convenient performance index has the quadratic form of 
J = j(x T Qx+w T Rw)dt ( 31 ) 



In order to design linear quadratic regulator the method of 
Lagrange multipliers, results in the following matrix Riccati 
equation 

PA + A T P + Q - PFR l F r P = (32) 
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Where, the pair of matrices (A F) and (A, C) are controllable 
and observable. 



The solution of (32) yields a positive definite symmetric 
matrix P and the desired optimal feedback gain matrix is given 
by 



K =R l F T P 



(33) 
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So the structure of optimal control law (30) is modifies as 

w = - R A F T Px (34) 
Using (28 and 30), A turns to A f 

A f = A- FK (35) 
And (28) turns to 



1 r ' 

J = -J"[(AC£,) 2 +(ACE 2 f +(ACE 3 f + (ACE 4 ) 2 + 
+ x 2 2 + x 22 2 ) + k(w 2 + w 2 
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+ w 3 2 + w 2 y\dt 



(37) 



(36) 

The Matlab control system toolbox function [K,P] = 
lqr2(A,F,Q,R) is used for the solution of matrix Riccati 
equation. 

The Performance Index selected for four-area power system is 
of the following form: 



If 2 



-a 34 x 16 + x ls + x 15 ) + (x 19 + x 20 + x 21 + x 22 ) 
+k(w 2 + w 2 + w 2 + w 2 y\dt 

From J, matrices Q and R can be recognized as given. 
Where, R (= KI) is Symmetric matrix. 
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Del Pd 1 = 0.01 



Del Pd 2 - 



Del Pd 3 - 



Del Pd 4 = 



x' = Ax+Bu 
y = Cx+Du 
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A Matlab program is written using matrices A, F, C, D, Q and 
R. The command[y, x] =step (Af, Fw, C, D, 1, t) is used to get 
the system's response. A Matlab Simulink model given in 
Fig. 3 can also be used to verify the system's response. 



Fig. 3 Simulink model of optimal four-area system using LQR design 

VI. RESULTS AND DISCUSSION 

Simulations were performed without controller, with 
Integral and the proposed optimal integral controllers applied 
to a four-area interconnected power system as shown in Fig.2, 
and Fig.3 by applying 0.01 puMW step load disturbance to 
Area 1. The same system parameters are used as given in [2], 
[11], [12]. The implementation is done on Matlab7.5 
Simulink- Workspace software. The simulations were run on a 
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personal computer Intel Core2Duo CPU T5450 @ 1.66 GH Z , 
982 MH Z , 2GB of RAM, under Window XP. The frequency 
deviations of Area 1 to Area 4 and some tie line power 
deviations with 0.01 puMW step load disturbance in Area 1 
are shown in Fig.4. The frequency deviation of Area 1 is 
presented in Fig.5 in a larger scale. The integral controller 



improves the system performance in comparison when no 
controller is used by eliminating steady state frequency errors. 
Moreover, the optimal integral controller is significantly 
superior to the integral controller. It gives a better performance 
than the integral controller. The settling time and peak 
overshoots are reduced considerably as shown in Fig.5. 



(a) 



(b) 

Fig. 4 Deviation in frequencies of (a) area 1, 2 and 3 (b) area 4, deviation in tie line power, AP 12 and AP lie , 41 for 1% change in load (APdi^O.01 pu). 
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Fig. 5 The frequency deviation of the first area in larger scale indicating peak overshoot and settling time for integral and optimal Integral controller for 1% 

change in load (APj^O.Ol pu). 
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The frequency deviation of Areal with settling time for 
5% band of the step load change and maximum overshoots is 
given in Table I. The comparison of dynamic performances of 
various controllers with the proposed optimal integral 
controller shows that the proposed controller gives better 
results in terms of lesser settling time and peak overshoot. The 
simulation was repeated with various instantaneous load 
changes and it was found that results from the controller are 
better every time. The simulation results show that proposed 
controller for load frequency control is given 47.68% 
reduction in settling time and 59.81% reduction in peak 
overshoots when compared with Cam and Kocaarslan study 
(for two-area with same data). 



TABLE I SYSTEM PERFORMANCE OF VARIOUS CONTROLLERS 



Deviation in frequency (Af i) 


Settling time (s) 


Peak over- 




for 5 % band of 


shoot (pu) 




step change 




Proposed study (Ta) 


1.8154 


-0.00779 


Cam's study (for two-area 


3.47 


-0.0194 


with FGPI) [11] 






Mathur's study ( for two-area 


3.53 


-0.0202 


with fuzzy controller) [10] 






Oysal's study (with NN) [15] 


6.81 


-0.0149 


Ghoshal's study (for three- 


8.0 


-0.028 


area with Fuzzy Controller) [9] 






Integral Controller (Tb) 


25.4803 


-0.0159 



In the analysis of the simulation results, the frequency 
deviation results were also used to calculate the integral 
absolute error (IAE) as given in Equation (39) for 40s of 
simulation time. 



IAE = J" I A/; I dt 

o 

TABLE II SYSTEM PERFORMANCE IN TERMS OF IAE 



Controllers 


IAE for Afi 


Optimal Integral Controller 


0.004663 


Integral Controller 


0.06361 


Uncontrolled System 


0.4757 



System performance of the controllers in terms of 
integral absolute error for frequency deviation of area 1 given 
in Table II also shows that proposed optimal integral 
controller has less integral absolute error as compared to 
conventional integral controller. 



VTI. CONCLUSION 

In this paper, an optimal integral controller is designed for 
automatic load frequency control of interconnected four-area 
power system by the use of optimal control design technique. 
The controller performance is observed on the basis of 
dynamic parameters i.e. settling time, peak overshoot and IAE. 
From Table I, it is concluded that the proposed controller for 
four-area has better dynamic performance when compared 
with many recent trend controllers applied to two-area and 
three-area interconnected power systems. 

APPENDIX 

The System data used in the study is given as: P r j= 2000 
MW, P rj = 2000 MW, T sgl = 0.08 sec, T psl = 20 sec, T tl = 0.03 
sec, R; = 2.4 H z pu/MW, T tJ = 0.08674, a§ = 1, K ti = 0.671, K psi 
= 120 H z pu/MW, b F 0.425 puMW/H z , AP dl =0.01 puMW, i, 
j= 1, 2, 3, & 4 



NOTATIONS 
th 

f ; nominal system frequency of i area, H z 

th 

Af; incremental frequency deviation of i area, H z pu 

th 

T sgi speed governor time constant of i area, sec. 

th 

R; speed governor regulation of the of i area, H z /puMW 

th 

T ti turbine time constant of i area, sec. 
Q positive semi-definite symmetric state cost weighting 
matrix 

R positive definite symmetric control cost weighting 
matrix 

k optimal feedback gain matrix 

th 

Kp Si gain of power system (generator load) of i area, 
H z /puMW 

th 

T psi power system time constant of i area, sec 

th 

APgi incremental generator power output change of i area, 
puMW 

th 

AP ti incremental turbine power output change of i area, 
puMW 

th 

AP ci incremental speed changer setting change of i area, 
puMW 

th 

AY Ei incremental steam valve setting change of i area, 
puMW 

th 

Kii gain of integral controller of i area 

th 

b; frequency bias of i area, puMW/H z 

th 

ACEj area control error of i area, puMW 

th 

Ui control input of i area 

th 

AP di incremental load demand change of i area, puMW 

th 

Wj disturbance vector of i area, puMW 
AP lie jj incremental tie line power change of i th and j th area, 
puMW 

ajj area size ratio coefficient 

Ty tie line synchronizing power co-efficient (pu) between 

area i and area j 
i & j subscript referred to area i & j. (i, j =1, 2, 3 & 4) 
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